The holographic induced gravity model with a Ricci dark energy: smoothing the little 

rip and big rip through Gauss-Bonnet effects? 
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We present an holographic brane-world model of the Dvali-Gabadadze-Porrati (DGP) scenario 
with and without a Gauss- Bonnet term (GB) in the bulk. We show that an holographic dark energy 
component with the Ricci scale as the infra-red cutoff can describe the late-time acceleration of the 
universe. In addition, we show that the dimensionless holographic parameter is very important in 
characterising the DGP branches, and in determining the behaviour of the Ricci dark energy as well 
as the asymptotic behaviour of the brane. On the one hand, in the DGP scenario the Ricci dark 
energy will exhibit a phantom-like behaviour with no big rip if the holographic parameter is strictly 
larger than 1/2. For smaller values, the brane hits a big rip or a little rip. On the other hand, we 
have shown that the introduction of the GB term avoids the big rip and little rip singularities on 
both branches but cannot avoid the appearance of a big freeze singularity for some values of the 
holographic parameter on the normal branch, however, these values are very unlikely because they 
lead to a very negative equation of state at the present and therefore we can speak in practice of 
singularity avoidance. At this regard, the equation of state parameter of the Ricci dark energy plays 
a crucial role, even more important than the GB parameter, in rejecting the parameter space where 
future singularities appear. 



I. INTRODUCTION 

Various astrophysical observations of type la super- 
novae yj , cosmic microwave background (CMB) and 
large scale structure (LSS) Q, etc., suggest that the uni- 
verse is in accelerating expansion. Quantitative analy- 
sis shows that there is a dark energy (DE) with negative 
pressure component leading to the dynamical mechanism 
of the accelerating expansion of the universe. However, 
the nature of this dark energy remains a mystery. Various 
models of dark energy have been proposed to solve this 
problem, such as a small positive cosmological constant 
4 1 and several kinds of scalar fields such as quintessence 
5|, k-essence Q, phantom etc. 

On the other hand, an alternative approach to explain 
the problem of DE arises from the holographic princi- 
ple which states that the number of degrees of freedom 
for a system within a finite region should be finite and 
bounded by the area of its boundary d, |^ . Furthermore 
it was also suggested that the ultraviolet (UV) cutoff 
scale of a system is connected to its infrared (IR) cut- 
off scale. Cohen et al. [l^] pointed out that, for a system 
with size L and UV cutoff, which is not a black hole, the 
quantum vacuum energy of the system should not exceed 
the mass of the same size black hole, i.e. L^pA < LMp, 
where p\ is the vacuum energy density caused by UV cut- 
off A, and Mp denotes the Plank mass (Af^ = I/SttG). 
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The largest IR cutoff L is required to saturate this in- 
equality [m 113 and its holographic energy density is 
given by pn = '^c^Mp/L^, where c is a phenomenologi- 
cal numerical constant. Later on, we will instead use the 
parameter j3 defined as /3 = . The holographic dark en- 
ergy model is based on assuming that the energy density 
Ph is responsible for the current speed up of the uni- 
verse with L being an appropriate cosmological length. 
It is well known that the length that plays the role of 
the IR cutoff on an holographic energy density is not 
unique. One of the choices of the IR cutoff is the Hubble 
rate. This choice does not induce acceleration either in a 
homogeneous and isotropic universe [l^ or in the DGP 
model [ll] , but it does in the DGP model with a Gauss- 
Bonnet (GB) term in the bulk [ll| (see also Ref. [Hi). In 
the reference , Li suggested that by choosing L as the 
event horizon, the late-time acceleration of the universe 
is suitably described. Another choice for the IR cutoff L 
was suggested by Gao ct al. [l^ (see also 17 1), in which 
the IR cutoff of the holographic Ricci dark energy (RDE) 
was taken to be the Ricci scalar curvature. 

Another approach to explain the observed acceleration 
of the late universe is to modify gravity at large scales. 
Models which describe this approach can be inspired by 
the existence of extra dimensions. The simplest and 
best studied case is the Dvali-Gabadadze-Porrati (DGP) 
mode l Il8i . Such a model contains two kinds of solu- 
tions [l9| : the self-accelerating branch, which suffers from 
some problems, and the normal branch. Even though the 
normal branch is healthy it cannot describe the current 
acceleration of the universe unless a dark energy compo- 
nent is invoked [l^ or the gravitational action is modified 
[2l| . Modifying gravity at large scales by introducing a 
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DGP brane-world model and taking into account curva- 
ture effects, as the GB term, were analysed in [22| , [11] , 
[23 | and [ij] where it was shown that they can be quite 
useful to describe the current acceleration of the universe. 
The main aim of this paper is to show if the holographic 
RDE is suitable to describe the late time acceleration of 
the universe, avoiding at the same time a big rip or little 
rip in a DGP model with a Gauss-Bonnet term in the 
bulk. 

The outline of the paper is as follows. In Sect. II 
we will review briefly the holographic DGP brane-world 
model with the Ricci scalar as an IR cutoif. In Sect. 
Ill, we analyse the modified Friedmann equation of the 
model without the Gauss-Bonnet term and we find an- 
alytically the asymptotical behaviour of the brane and 
numerically the whole expansion of the brane. In par- 
ticular, we discuss the behaviour of the equation of state 
associated with the holographic Ricci dark energy and 
we show that for /? < 1/2, /3 being the holographic pa- 
rameter, the solutions lead to a big rip or little rip. This 
motivate us to consider a GB term in the bulk action in 
order to enlarge the /3 spectrum and overcome this prob- 
lem. In Sect. IV, we consider the model where the bulk 
contains a GB term. We show that the effect of a GB 
term is to remove the little rip and big rip singularities 
on both branches, but it cannot avoid the appearance of 
a big freeze singularity for some values of the holographic 
parameter f3 on the normal branch, however, these val- 
ues are very unlikely because they lead to a very nega- 
tive equation of state at the present and therefore we can 
speak in practice of singularity avoidance. At this regard, 
the equation of state parameter of the Ricci dark energy 
plays a crucial role, even more important than the GB 
parameter, in rejecting the parameter space where future 
singularities appear. Finally, in Sect. V, we conclude. 



the GB parameter, respectively, both of them being pos- 
itive. The parameter e in Eq. (|2.ip can take two values: 
e ~ 1, corresponding to the self- accelerating branch in the 
absence of any kind of dark energy [2^ ; and e = — 1 , cor- 
responding to the normal branch which requires a dark 
energy component to accelerate at late-time [l3, HH, [13] . 
For simplicity, we will keep the terminology: (i) self- 
accelerating branch when e = 1 and (ii) normal branch 
when e = — 1. 

As already mentioned in the introduction, is related 
to the UV cutoff, while L is related to the IR cutoff. 
We consider a Ricci dark energy (RDE) [l^ where the 
scale L is fixed by TZ/6, TZ being the scalar curvature 
of the brane. For a spatially flat Friedmann-Lemaitre- 
Robcrtson- Walker (FLRW) universe 



n=6{H + 2H' 



where the dot stands for the derivative with respect to 
the cosmic time of the brane. Therefore, the RDE energy 
density is 



PH = 3/3 A/^ - 



1 dH^ 



2 dx 



2W 



(2.2) 



where x = — ln(z+l) = ln(a), z is the rcdshift and j3 = 
is a dimcnsionless parameter that measures the strength 
of the holographic component. As we will show, /? and 
Eq. (j2.ip play a crucial role in determining the expan- 
sion of the RDE component. Before moving forwards, 
we would like to highlight that we have assumed that the 
holographic dark energy stated in Eq. (|2.2|) remains valid 
on brane-world models like the one we are considering 
here. A similar approach was followed in Refs. [l3l. IT^. 

The modified Friedmann equation (|2.ip can be further 
rewritten as: 



II. 



SETUP OF THE MODEL 



We consider a DGP brane-world model, where the bulk 
action contains a GB curvature term. The bulk corre- 
sponds to two symmetric pieces of a 5-dimensional (5d) 
Minkowski space-time. The brane is spatially flat and its 
action contains an induced gravity term. We assume that 
the brane is filled with matter and an holographic dark 
energy. Then, the modified Friedmann equation reads 
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(2.1) 



1 dE'^ 

= 17™e-3^ + /3( + 2E^) + 2eJnZ{l + ^aE^)E, 

2 dx 



where E{z) = H/Hq and 

Pmo 



(2.3) 
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4^' and = -aH', 



(2.4) 



are the usual convenient dimensionless parameters and 
the subscripts denotes the present value (we will follow 
the same notation as in [3: [13 ) ■ The dimensionless pa- 
rameter of the holographic energy density can be written 
as: 



where H is the brane Hubble parameter, p ~ p,n + pn is 
the total cosmic fluid energy density of the brane which 
can be described through a cold dark matter component 
(CDM) with energy density p„i and an holographic Ricci 
dark energy component with energy density pu • The pa- 
rameters Tc and a correspond to the cross over scale and 
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From (|2.3p the variation of the dimensionless Hubble rate 
E with respect to x is 
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dE 
dx 



nr. 
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(2/3 -i)E^ + 2eVn;:;(i + naE^)E 



I3E 



(2.6) 



It is convenient to rewrite the previous equation as 



E 



+ (2/3-1)^2 



(2.7) 

From this equation, we can conclude that: (i) on the 
one hand, H is well defined and has no singularity at 
any finite dimensionless Hubble rate. This indicates the 
absence of sudden singularities f26l - [29| in this model, (ii) 
On the other hand, if E diverges, the energy density pa 
(cf. Eq. (|2.2p ) could also diverge and the brane could 
undergo a big rip or a big freeze singularity [H, [s^, HH 
or even a little rip [32| - [33 | . This issue will become clearer 
in the next sections. 

Even though Eq. (|2.7p describes the first derivative of 
the Hubble rate, it is not the Raychaudhuri equation but 
the Friedmann equation (|2.3p . The appearance of H in 
Eqs. (|2.3|) or (|2.7p is simply a consequence of the defini- 
tion of the holographic dark energy (cf. Eq. (|2.2|) ) which 
depends on the first derivative of the Hubble rate. The 
Raychaudhuri equation can be easily obtained from the 
modified Friedmann equation (|2.3|) and the conservation 
of the matter content of the brane. We skip it here as 
the equation is very lengthly and not very enlightening. 

The modified Friedmann equation ()2.3|) evaluated at 
the present time implies a constraint on the set of pa- 
rameters of the model 
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In addition, we have a couple of constraints on the dif- 
ferential equation (|2.6p 



Eo = E{x 
dE/dx\^^o = 



0) 

-2- 



= 1, 



(2.9) 



where the second condition comes from the dimension- 
less energy density definition (see Eq. (|2.5I) ) evaluated at 
the present time. It is possible to constrain further the 
cosmological parameter by using the deceleration param- 
eter, q= — (1 + ^dE/dx), which reads at present: 



90 = 1 - 



(2.10) 



The universe is currently accelerating, i.e. go < 0, and 
the holographic RDE energy density is positive, conse- 
quently 

0</3<f2Ho- (2.11) 

An analytical solution of Eqs. (|2.6p or (|2.7p is not at 
all obvious, so we will combine an analytical asymptoti- 
cal analysis with a numerical one to study the modified 
Friedmann equation (|2.6p . We will also split our analysis 
in two parts: (i) a purely induced gravity brane- world 
model and (ii) the same kind of model with GB correc- 
tions in the bulk. 
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FIG. 1. Plot of the dimensionless Hubble rate E against the 
redshift 2 for the normal branch (0.43 < P) and for the self- 
accelerating branch (/3 < 0.43). The cosmological parameters 
are assumed to be ilm ~ 0.27 and go = —0.7. The plot has 
been done for several values of the holographic parameter /? 
that are stated at the bottom of the figure. 



III. HOLOGRAPHIC RICCI DARK ENERGY IN 
A DGP BRANE- WORLD 

In order to continue the analysis of this model it is con- 
venient to start with the simplest case for which a = 0. 

As we know, the universe is currently accelerating, i.e. 
qa < 0. Hence, using Eqs. ((2^ and ((279| (for a = 0), we 
obtain 



</3 



1-a. 



(3.1) 



This inequality is quite enlightening and implies an upper 
or a lower bound on the allowed values of /3 depending 
on which DGP branch we are considering: 



where 



13 < Aim 
/? > Aim 



/3lim — 



for 
for 



e = -l, 



1 i^rn 

1 - go 



(3.2) 



(3.3) 



We can estimate such a bound. Aim, taking into account 
the latest WMAP7 data Q and considering that our 
model does not deviate too much from the ACDM sce- 
nario at the present time; i.e. ilm ~ 0.27 and go ~ —0.7, 
and so we get a threshold value Aim ~ 0.43. From now 
on, for a given value of /?, we know with which DGP 
branch we are dealing with, depending on whether /S is 
smaller or larger than 0.43. 

As wc mentioned at the end of Sect. [TTl it is not at all 
obvious how to solve analytically the Friedmann equa- 
tion (|2.6p . even in the simplest case {a = 0). We will. 
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FIG. 2. Plot of the equation of state parameter u-^ against the 
redshift z for the normal branch (0.43 < P) and for the self- 
accelerating branch (/3 < 0.43). The cosmological parameters 
are assumed to be Q,m = 0.27 and go = —0.7. The plot has 
been done for several values of the holographic parameter /? 
that are stated at the bottom of the figure. 



FIG. 3. Plot of the deceleration parameter q against the red- 
shift z for the normal branch (0.43 < /3) and for the self- 
accelerating branch (/3 < 0.43). The cosmological parameters 
are assumed to be Q,m = 0.27 and go = —0.7. The plot has 
been done for several values of the holographic parameter /? 
that are stated at the bottom of the figure. 



therefore, analyse the asymptotical solutions of Eq. (|2.6|) 
and combine them with numerical solutions of the same 
equation. In addition, this analysis is very important to 
analyse the asymptotical behaviour of the brane and the 
possible occurrence of any doomsday in the future. 

In the far future (z — >■ —1), the energy density of mat- 
ter = Pmai^ + zY is practically zero, and the universe 
converges asymptotically to a universe filled exclusively 
with an holographic RDE component. 

In the absence of matter, we distinguish two cases for 
the solutions of equation (|2.6p with a = 0: 



A. Asymptotic behaviour (3 = 1/2 

We restrict our analysis to expanding branes. There- 
fore, there is a unique solution supported by the normal 
branch as implied by the inequalities p.2p : 

^ (3.4) 



where Ci is an integration constant. The solution p. 41 
can be further integrated over time 



1 



4yn^-ffo(i - ti) 



exp 



Ci}, (3.5) 



where xi and ti are integration constants. Consequently, 
in the far future, the scale factor (or x) and the Hubble 
rate diverge. What about HI It can be easily shown 
that 



dE 
'dt 



4 a/17, 



(3.6) 



and therefore the cosmic derivative of the Hubble rate 
also diverges. This kind of event has been named recently 
a little rip f33 | and was previously found on a 4D model 
(3^ and also on a dilatonic brane- world model [s^. In 
the latter case, the event was called a big rip singularity 
sent towards an infinite future cosmic time 3^ . 



B. Asymptotic behaviour /3 7^ 1/2 

The solutions can be expressed as: 




exp 



2(3- 1 



{x - Xi) 



(3.7) 



where Ei and Xi are integration constants, and can be 
split in two subsets: 



1. 1/2 < p 

This set of values for (3 corresponds to the normal 
branch as discussed just after Eq. (|3.2|) . and the dimen- 
sionless Hubble rate takes the solution 



E+ = 



(3.8) 



(2/3-1) 

and becomes asymptotically de Sitter in the far future. 
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2. 



P<l/2 



This case can correspond to any of the two DGP 
branches (please see the argument given just after 
Eq. p.2p ). The dimensionless Hubble rate blows up in 
the far future at a finite future cosmic time, and it fol- 
lows a super-accelerated expansion until it hits a big rip 
singularity [28l. l30l|. 



C. Numerical analysis 



Now we return to the equation (j2.6p with a = 0, in 
order to solve it and to study the behaviour of the holo- 
graphic Ricci dark energy. As already mentioned, this 
equation has no analytical solutions, and a numerical 
analysis is required. We show our results in Figs. [U [2] 
and[3l 

In order to solve numerically Eq. (|2.6p . we choose 
~ 0.27 and qo ^ -0.7, as given by the latest WMAP7 
data and assuming that our model is very similar to 
the ACDM scenario at the present time. For a given 
value of /?, the dimensionless energy density fiuo is fixed 
through Eq. (|2.10p . while the crossover scale parameter 
fi^c is fixed by the constraint equation (|2.8p . 

In Fig. [U we show five solutions of Eq. (|2.6p corre- 
sponding to different values of /3. All the solutions start 
with a big bang singularity while their behaviour in the 
future depends strongly on the values acquired by j3 (cf. 
Fig.[T]): for /3 = 0.4 and /3 = 0.45, the brane hits a big rip 
singularity in the future, for j3 = 0.5 the brane expands 
during an infinite time until it reaches a little rip, and 
finally for /3 = 0.6 and /3 = 0.8 the brane is asymptoti- 
cally de Sitter in the future and the dimensionless Hub- 
ble rate E approaches the value = ^^^~(2/3-Ii)^""' (^^^ 
Eq. p.8|) ). All these numerical solutions are in agreement 
with our analytical analysis as should be. 

Even though we have already imposed that the brane 
is currently accelerating, more precisely we imposed 
go —0.7, it is important to analyse the behaviour of the 
equation of state parameter, wh = Ph/ph, and the decel- 
eration parameter, q, as functions of the redshift. Using 
the equation of conservation of the Ricci holographic dark 
energy (HDE), 



Ph + 3i?(pH + Ph) = 0, 



we obtain 



-1 



1 din Oh 

3 dx 



(3.9) 



(3.10) 



r^H is defined in equation (|2.5p . In terms of the dimen- 
sionless Hubble rate E and its derivatives with respect 
to X, wh can be rewritten as: 

2 



-1 



dE 
dx 



E' 



J'E 
dx"^ 



AE 



dE 
dx 



+ 6E^ 

dx 



(3.11) 



In Fig. [2l we show an example of the behaviour of the 
equation of state parameter, ojh, for the current cosmo- 
logical values, and for different values of (3. 

At very late-time wh approaches —1 for /? > 1/2 and 
tends to ^^-^ for l3 < 1/2 as can be clearly proved from 
Eq. ((^ and Eq. ([XT)) . We would like to highlight in 
this regard the smooth transition on the values of /3 when 
going from an asymptotically de Sitter brane (1/2 < f3) 
to a little rip (/? = 1/2) to finally a big rip {13 < 1/2). 
In addition, in this model we can see the little rip as 
a division between an asymptotically dc Sitter universe 
and a universe approaching/hitting a big rip singularity. 

On the other hand, one can prove that for all values of 
/3 the equation of state of dark energy will evolve in the 
future in the region cjh < ^1- In addition, using equa- 
tions (|2.6p and (|3.1ip , the parameter wh can be expressed 
at the present time as: 



\(3 (qg - 1) - (go-2)-(l-t-go) 
^Ho = ^^jt; ^ • (3.12) 



3/3(1 -go) 



From the current observational data and from (|2.1ip 
we can assume that /3 < 0.73, which means that < 
— 1, and the Ricci dark energy has a phantom-like be- 
haviour. For /3 > 0.5, this behaviour will continue in 
the future and will be more and more like a cosmological 
constant such that ultimately the universe will enter a 
de Sitter phase in the far future. For /3 < 0.5 the be- 
havior of the Ricci HDE will be phantom-like until the 
brane reaches a big rip or a little rip. This behaviour is 
not consistent with the holographic principle which pro- 
hibits the occurrence of a big rip singularity (25l |. This 
motivates us, as was already done in the case of an holo- 
graphic dark energy with UV cutoff L ^ [l3l, to 
introduce a GB term in the bulk to try to remove this 
singularity. 

Similarly we can plot the deceleration parameter q 
against the redshift z for different values of /3, and for the 
current cosmological values. The curves in Fig. [3] show 
that our universe has entered a phase of accelerated ex- 
pansion in the recent past. The deceleration parameter 
reaches at very late-time the value —1 for 1/2 < /3, and 
the value ^ for ;3 < 1/2. 



IV. HOLOGRAPHIC RICCI DARK ENERGY IN 
A DGP BRANE- WORLD WITH A 
GAUSS-BONNET TERM IN THE BULK 

We consider now the model where the bulk action con- 
tains a Gauss-Bonnet (GB) curvature term in order to 
analyse the possibility of avoiding the little rip and big 
rip present in the absence of such a term. Our motivation 
in including this curvature modification is based on the 
fact that a GB curvature term corresponds to an ultra- 
violet correction to gravity and therefore can influence 
the high energy regime where both the little rip and big 
rip take place. 
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FIG. 4. Plot of the parameters /3± and /3iim, defined in 
Eqs. (|4.5p and (|3.3p . respectively, versus the parameter JIq. 
We have used the current observational data 50 ~ —0.7, and 
rJm 0.27. The parameter /Jum defines the border line be- 
tween the normal branch (/3iim < P) and the self-accelerating 
branch (/3 < /3iim). For a given value of /3, the parameters /?± 
fix the sign of the discriminant T) (cf Eq. (|4.4p 'l. 



By combining the expressions 
tain: 



and (|2.10p . wc ob- 



1 - n„, - - go) 
2(1 + 17a) 



(4.1) 



Therefore, the right hand side (rhs) of the previous equal- 
ity is positive for the self-accelerating branch (e = 1) and 
negative for the normal branch (e = —1). Consequently, 
there is a limiting allowed value for /3, given by /3iim and 
defined in Eq. ()3.3|) . depending on which branch wc arc 
considering such that the inequality p.2p still holds. Fol- 
lowing the same procedure introduced in Sect. IIIIl we 
conclude that /3iim ^ 0.43. In summary, a value of /3 
such that /3 < 0.43 corresponds to the self-accelerating 
branch, while a value of /3 such that 0.43 < /3 corresponds 
to the normal branch. 

We would like to highlight that this model has two 
free parameters: for example fie and /3. Indeed, for a 
given set of values (12q,/3), the model is fully defined be- 
cause is fixed by observations, r^Ho is determined by 
Eq. (|2.10p and $7^^ by Eq. (|4.ip . In our reasoning we have 
assumed that the ACDM scenario is a good approxima- 
tion of our model at the present time and we have used 
the latest WMAP7 data 0. 

With the inclusion of a GB term in the bulk, it is 
even less obvious how to solve analytically the Friedmann 
equation (|2.3p . We will, therefore, analyse the asymptot- 
ical solutions of Eq. (|2.3p (or equivalently Eq. (|2.6p ) and 
combine them with numerical solutions of the same equa- 
tion. This analysis will clarify the possible avoidance of 
a doomsday, of the kind of big rip or little rip, present on 
the same scenario without GB effects for a set of values 
of /3. 

In the far future (z — >■ —1), the energy density of mat- 
ter pm ~ Prnoi^ + -z)"^ cau be neglected, and the universe 
converges asymptotically to a universe filled exclusively 
with an holographic RDE component. 



A. Asymptotic behaviour 

In the absence of matter, the modified Friedmann 
equation (|2.6p can be rewritten as: 



/3 



dE 

dx 



{213 ~ I) E + 2e^/nZ{l + ^c^E'^) , (4.2) 



whose solution depends on the sign of the discriminant 
T) of the polynomial on the rhs, which reads 

V = {213 ~lf -IQVLrna. (4.3) 
It is convenient to factorise the discriminant T) as follows 



V = T{P-P^){P-p+). 



where 



/3± 



1 -h f7„ ± 2\/n^(i - n,„) 

2 [1 + ± vni(i - go)] ' 
2" 



J" = 4 



1 - Via 



1 - go 

l + f7a 



(4.4) 

(4.5) 
(4.6) 



We show in Fig. 31 the shape of (3± for physically reason- 
able values of Q-a- It can be shown as well that for the 
same range of values of Via, the coefhcient F is positive 
(see Eq. (|4.6p ). Therefore, the discriminant T) is positive 
when /3 < /3+ or /3_ < /3, negative for (3+ < 13 < j3- and 
vanishes for /3 = , /3_ . We next tackle each of these 
cases separately. 



1. Negative discriminant (D <Q) 

The discriminant T) is negative when (3 satisfies /?+ < 
/3 < /?_ (cf. Eq. (gH) and Fig. S]). As Aim < 13+ this 
case can take place only for the normal branch and the 
Friedmann equation (|4.2p can be integrated as: 




and X2 and E2 are integration constants. Consequently, 
there is always a finite value of the scale factor or x where 
the Hubble rate blows up: 



''singi 




2B 1 

(n+-)7r, 71 eZ. (4.9) 



The singularity takes place at Xsi^g for an integer n such 
that Xsingi takes the smallest possible positive value. Not 
only the Hubble rate but also its cosmic derivative blow 
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where 



1 - \/P 



4en„yn;~£'2 + 2/3 - 1 + \/P 



(4.12) 



X2 and El are integration constants. We can then deduce 
that the brane expands till it reaches an infinite size; i.e. 
for very large values of a;, if and only if the brane is 
asymptotically dc Sitter with Hubble rate 



V-{2f3- 1) 



(4.13) 



I p-0.45 p-0.7 |i-O.S P-0.4| 



FIG. 5. Plot of the dimensionless Hubble rate E against x 
for the normal branch (0.43 < /3) and for the self-accelerating 
branch (/? < 0.43). The cosmological parameters are assumed 
to be Qrn = 0.27 and qo = —0.7. The plot has been done for 
several values of the holographic parameter /3 that are stated 
at the bottom of the figure. The reason we show this plot in 
terms of x, rather than z, is to see clearly the presence of a 
big freeze singularity for some values of /3. 



up at Xsingi (see Eq. (|2.7|) '). Finally, integrating the di- 
mensionless Hubble rate (|4.7p with respect to time we 
obtain 



X>|ln< 



y|p[tan 




+ 2/3-1 


\l 


1 + tan^ 


2/3 ^2) 





Ho + (2/3 - 1)2] 



(i ^singi); 



> + 2/? - 1 



(4.10) 



Notice that i?+ reduces to Eq. (|3.8p in absence of GB 
effects in the bulk. The brane can be asymptotically de 
Sitter only in too cases: (i) For e = 1 as /3 < /3ii,n < \ 
and (ii) for e = — 1 if ^ < /3_ < /3. In both cases the 
Hubble rate £'+ is a positive constant. 

What happens in the remaining case; i.e. e = — 1 and 
Aim < /3 < /3-I-? In the latter case Ej^ is negative and the 
asymptotic analysis performed is no longer valid. In fact, 
the brane faces a big freeze singularity in this situation 
where the expansion of the brane can be approximated 
by 



E 



1 



X) 



(4.14) 



The constant a;sing2 stands for the "size" of the brane at 
the big freeze singularity. The Hubble rate (|4.14p can be 
integrated over time, resulting in the following expansion 
for the scale factor of the brane 



exp 



sing2 



- t 



(4.15) 



The big freeze singularity takes place at a finite scale 



factor, 



and a finite cosmic time, 



concluding therefore that the brane hits a big freeze sin- 
gularity in the future as the event Xsingi takes place at a 
finite future cosmic time igingi [H, |3l| . 



2,. Positive discriminant (0 <T>) 

The discriminant V is positive when /3 < /3_|_ or /3_ < /3 
(cf. Eq. (j4.4p and Fig. |4]) which can take place either on 
the normal branch (/3 > /3iim) or the self-accelerating 
branch (/3iini > /3). 

The solution of the asymptotic Friedmann equation 
(14:21) reads: 



1 - 



AenaJTl^^E + 2/3 - 1 + 



= C2 cxp 



- — (.X-X2) 



3. Vanishing discriminant (D = 0) 

Finally, the discriminant V vanishes when /3 = /3_)_ or 
/3_. This can take place only on the normal branch. In 
this fine tuned case, the solution of the modified Fried- 
mann equation (j4.2p can be expressed as 



E 



where 



(2/3-1) 



1 



C2 + 2^(x - X2) 



(4.16) 



(4.17) 



and X2 and E2 are integration constants. By taking the 
limit a; — > 00, the dimensionless Hubble rate (|4.16p re- 
duces to 



E^ = 



2/3-1 



(4.18) 



\ 

I 
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FIG. 6. Parameter cu^ against x for the normal branch 
(0.43 < P) and for the self-accelerating branch (/3 < 0.43). 
The cosmological parameters are assumed to be = 0.27, 
go = —0.7 and Qa ~ 0.1. The plot has been done for several 
values of the holographic parameter /3 that are stated at the 
bottom of the figure. The reason we show this plot in terms 
of X, rather than z, is to see what are the effects of the big 
freeze on the equation of state parameter. 



Not surprisingly, coincides with the one defined in 
Eq. (|4.13p for vanishing V on the normal branch. The 
dimcnsionless Hubble rate is positive only for values of (3 
such that /3 = /?_ where the brane is asymptotically de 
Sitter. The solution (j4.16p is not valid on the normal ex- 
panding branch with /? = , as the Hubble rate becomes 
negative, a condition which signals that the asymptoti- 
cal behaviour given by (|4.16p breaks down. In fact, what 
happens on the normal branch for /3 = /?_|_ is that the 
brane hits a big freeze singularity at a finite scale fac- 
tor and the brane expansion can be approximated by 

Eq. mn) . 



B. Numerical analysis 

The asymptotical analysis we have carried out in the 
previous subsection can be complemented with a numer- 
ical analysis. We show some examples of the numerical 
analysis in Fig [5] More precisely, we obtain: 

• If /3 < Prnn; the solution corresponds to the self- 
accelerating branch and the brane is asymptotically 
de Sitter in the future. 

• If Piim </?</?-, the solution corresponds to the 
normal branch and the brane faces a big freeze in 
the future. Notice that for the case /3 = 1/2 with 
r^Q, 7^ the brane faces a big freeze, while it faces 
a little rip in the absence of a GB term in the bulk 

(f^a = 0). 



FIG. 7. Plot of iJHo against /3 and Qa ■ The cosmological 
parameters are assumed to be Q,m ~ 0.27 and go = —0.7. 

• If /3_ < 13, the solution corresponds to the normal 
branch and is asymptotically de Sitter in the future. 

In summary, we conclude that the big rip that might 
show up on the self-accelerating branch without GB ef- 
fects in the bulk is removed by the inclusion of GB effects 
in the bulk. However, the situation in the normal branch 
is quite different as the presence of the GB effects, even 
though they remove the big rip and little rip, it cannot 
always prevent the presence of a big freeze singularity in 
the future evolution of the brane. 

For completeness, we will discuss the behaviour of the 
equation of state wh associated with the Ricci HDE. We 
can show that the function wh approaches —1 at very 
late-time, in the case where the brane is asymptotically 
de Sitter, while it diverges to minus infinity when the 
brane face a big freeze singularity in the future (see also 
Fig.ini). In addition, we can show by using equations (|2.6p 
and p. lip , that the equation of state parameter wh can 
be expressed at the present time as: 

^Ho = 

[/3{1 - go) + [go(l + Sn^) - 2] - (1 + go) {Qa - 1) 
3/3(1- go) (f^a + 1) 

(4.19) 

Fig. [7] shows the variation of ojho with respect to fla and 
/3, where we have used the current values of go ^ —0.7 
and ~ 0.27. The figure shows that is always less 
than —1/3 at present. Notice as well that the most likely 
value of (3 is of the order of (3- or larger, and therefore 
the big freeze singularity will be avoided in the future 
evolution of the brane (cf. Figs. |4] and [7]). 

One can conclude that for all values of fta and /3, the 
equation of state of dark energy wh will evolve at future 
in the region cjh < —1/3, and converges to —1 or minus 
infinity as discussed above. Fig. IH] illustrates the case 
of ilct = 0.1. This figure shows that the curve of wh 
crosses currently the boundary of the cosmological con- 
stant wh = — 1 for /3o ~ 0.77, while for (3 < /3o the equa- 
tion of state of dark energy will evolve as a phantom-like 
fluid, i.e. wh < — 1. This behaviour will manifest itself 
more and more in the future until ending with a big freeze 
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FIG. 8. Plot of the deceleration parameter against x for the 
normal branch (0.43 < /3) and for the self-accelerating branch 
(/J < 0.43). The cosmological parameters are assumed to be 
= 0.27, go = -0.7 and = 0.1. The plot has been 
done for several values of the holographic parameter /? that 
are stated at the bottom of the figure. The reason we show 
this plot in terms of x, rather than 2:, is to see what are the 
effects of the big freeze on the decceleration parameter. 



for /3iim</3 < /?-, and with a dc Sittcr-like universe for 
/3 < Aim and /?_ < /?. 

We note here that according to Eq. (|2.11[) /3 is always 
less than /3o and the nature of the Ricci dark energy has 
always a phantom-like behaviour. 

Once again one can test this model by plotting the de- 
celeration parameter q versus the redshift for some values 
of /3. Fig. [5] shows that our universe is currently in an 
accelerating expansion phase. 



V. CONCLUSIONS 

In this paper we have considered an holographic DGP 
brane-world model with the average radius of the Ricci 
scalar curvature as an infrared cutoff. We have shown 
that the late-time acceleration of the Universe with and 
without the Gauss-Bonnet (GB) term in the bulk, is con- 
sistent with the current observational data. 

The Ricci HDE is characterised by the parameter /3 
(see Eq. (|2.2p ) which is bounded by the current value of 
the dimensionless Ricci HDE (cf. Eq. (|2.11|) ). In 
addition, there is a limiting value of /?, which we named 
Aim (see Eqs. (|3.2p and (|3.3p ) that splits the allowed /3 
values for the self-accelerating branch (e = 1) from the 
normal branch (e = —1). More precisely, values of /3 such 
that /3 < Aim correspond to the self-accelerating branch 
while values of /3 such that Aim < /3 corresponds to the 
normal branch. We can estimate that Aim ^0.43 assum- 
ing that our model does not deviate too much from the 
ACDM model at present. The value acquired by /3 is very 
important in determining the asymptotical behaviour of 



the brane and that of the Ricci HDE. 

In the absence of GB effects, the equation of state of 
the Ricci HDE will evolve in the region wh < — 1 and 
its late-time behaviour is phantom-like with two distinct 
destinies of the universe: (i) for 1/2 < /3, correspond- 
ing exclusively to the normal branch, this phantom-like 
behaviour will be more and more like a cosmological con- 
stant in the far future and the brane expansion is asymp- 
totically dc Sitter. In this case, the normal branch docs 
not require the inclusion of a GB term in the bulk to 
remove any singularity in the far future, (ii) for /? < 1/2 
(self-accelerating or normal branch), this phantom-like 
behaviour continues and implies a super-accelerated ex- 
pansion on the brane leading to a big rip or a little rip 
singularity. These sets of solutions may require a GB 
term in order to remove or smooth those singularities. 
This is what motivates us to consider a GB correction in 
the bulk. 

The introduction of a GB term in the bulk may re- 
move the big rip and little rip from the future evolution 
of the brane. In the self-accelerating branch the brane 
becomes asymptotically de Sitter, however, the normal 
branch can hit a big freeze singularity in the far future. 
More precisely, the normal branch is asymptotically de 
Sitter for /3_ < /3 and hits a big freeze singularity when 
Aim < P < 13- (see Eqs. (gS]), (4.7) and (gill)). The 
Ricci HDE has a phantom-like behaviour which either 
approaches a de Sitter-like state or induces a big freeze as 
explained above. The case corresponding to a big freeze 
singularity is quite unlikely because it leads to a quite 
negative equation of state at present (please see Figs. |3] 
and[7l). At this regard, we would like to stress that an 
observationally consistent value of the equation of state 
(EOS) parameter of the Ricci HDE is even more impor- 
tant than the GB parameter to reject the parameter re- 
gion where the future singularity appears. In fact, even 
in the absence of a GB parameter, the EOS parameter al- 
ready pick up a singularity free region (cf. Fig. 2). Then 
what is the role of the GB term? it narrows the region af- 
fected by the singularity if /3 < 1/2 (with respect to the 
model without a GB term) and enlarges that region if 
1/2 < /3 (with respect to the model without a GB term). 

In summary, it is always possible to construct a Ricci 
HDE brane-world model consistent with current observa- 
tions. The brane with or without GB effects is always free 
of future singularities except for the range of values of (3 
such that Aim < /3 < 1/2 which happens always in the 
normal branch. However, these values are very unlikely 
(cf. Figs. |4] and ^ because they lead to a very negative 
equation of state at the present as explained above. 

This study opens a new way of considering interactions 
between dark matter and the holographic dark energy. 
We hope to report this issue in a forthcoming publication 
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